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ABSTRACT 


An  experimental  natural  convection  heat  transfer  study  of  a  simulated  electronic 
circuit  board  has  been  conducted.  The  board  has  an  array  of  9  simulated  chips  ,  each 
dissipating  up  to  2.5  Watts.  The  board  is  immersed  in  FC75.  a  fluorocarbon  liquid,  in 
an  enclosure  whose  top  and  bottom  surfaces  are  constant  temperature  heat  sinks.  The 
experimental  data  have  been  expressed  in  terms  of  relevant  dimensionless  heat  transfer 
parameters  such  as  Nusselt  and  Rayleigh  numbers.  The  trend  is  that,  the  chips  located 
higher  in  the  enclosure  have  lower  heat  transfer  rates.  Otherwise  the  chips  in  the  same 
row  behave  in  a  similar  way  which  implies  a  quasi-two  dimensionality. 
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I.  INTRODUCTION 


A.  DESCRIPTION  OF  THE  PROBLEM 

Development  of  electronic  packaging,  starting  with  ENIAC,  the  first  digital 
electronic  computer,  has  created  a  very  urgent  need  to  devise  more  sophisticated 
cooling  systems.  In  the  last  ten  years,  the  remarkable  advancement  of  very  densely 
packaged  electronic  components  has  brought  about  the  necessity  for  enhancement  of 
heat  transfer  features  applied  to  these  thermal  packages. 

The  first  electronic  packages  in  I960's,  SSI  (Small  Scale  Integration)  chips  had 
only  a  maximum  of  20  devices  per  chip,  whereas  today's  most  advanced  L'LSI  (Ultra 
High  Scale  Integration)  chips  have  more  than  100,000  devices  per  chip  [Ref.  1.]  The 
overall  effect  over  the  past  30  years  resulting  from  this  tremendous  increase  in  the 
packaging  densities  has  been  a  trend  towards  ever  increasing  heat  fluxes.  While  this  has 
always  been  the  case  for  very  high  power  modules,  the  same  trend  towards  very  high 
heat  fluxes  in  the  microelectronic  chips  has  also  been  observed.  Heat  fluxes  reaching 
105  W  m2  on  the  chip  surface  have  introduced  a  major  cooling  problem  [Ref.  2J. 
Inadequate  cooling  of  these  elements,  coupled  with  the  enviromental  temperature  rise, 
often  leads  to  failure  of  these  components  and  decrease  in  reliability  and  performance. 

B.  BACKGROUND 

Many  papers  on  the  electronic  cooling  area  have  proposed  a  very’  wide  variety  of 
solutions  to  specific  problems.  Major  digital  technology  firms  such  as  IBM,  Honeywell 
and  Mitsubishi  have  become  the  pioneers  in  this  research  area.  For  instance,  air 
cooling  technology  of  the  modules  now  often  includes  impingement  air  cooling  to 
enhance  the  heat  transfer.  Using  liquid  cooled  cold  plates  to  conduct  heat  away  from 
the  sources  has  also  been  employed.  Combined  conduction  and  air  cooling  has  been 
observed  to  provide  a  much  higher  enhancement  when  used  with  fins  on  the  chip 
surface  [Ref.  3]. 

Some  special  chip  outer  surface  features  such  as  TPF  (Turbulence  Promoting 
Fin)  have  provided  significantly  higher  cooling  enhancements,  when  used  with  special 
gases  such  as  Helium  as  the  closed  cycle  coolant  [Ref.  4],  By  doing  so.  300%  greater 
dissipation  (cooling)  level  can  be  reached  than  when  open  cycle  air  is  used  as  the 
coolant. 


tMNUlOMNUV'UV 


a 

1 


U 


iWWVWW, 


Among  the  cooling  techniques  mentioned  above,  liquid  forced  convection  has  a 
very  important  place  since  it  is  possible  to  attain  heat  transfer  coefficients  over  an 
order  of  magnitude  higher  than  those  attainable  with  air  cooling  [Ref.  1J.  According  to 
TIME  magazine.  June  17,  19S5  issue,  the  Cray-2  super  computer  utilizes  forced 
convection  with  local  boiling  at  some  hot  spots.  In  the  same  magazine  it  is  stressed 
that  "2u0  gallons  of  liquid  arc  used  in  the  machine  flooding  the  circuits  with  a 
continuous  flow  of  Fluorinert."  However,  the  use  of  forced  convection  leads  to  more 
complex  design  features  and  brings  along  the  possibility  of  vibration  based  failure  risks. 

Finally,  more  recent  advancements  and  growing  use  of  largescalc  and  very  large 
scale  integration  (LSI  and  VLSI)  microelectronic  technologies  with  chip  dissipation 
exceeding  1-2  VV  and  chip  surface  heat  fluxes  exceeding  10s  W  m2  has  focused 
attention  on  natural  convection  liquid  immersion  cooling  because  of  its  advantages 
over  forced  convection  such  as  simplicity,  higher  reliability  and  lower  cost  [Ref.  5]. 
Using  R- 113  as  the  coolant,  different  immersion  cooling  schemes  with  submerged 
condensers  have  been  proposed  by  3ar-Cohen  [Ref.  6].  These  techniques  are  mainiy 
based  on  pool  boiling  of  R- 113  or  other  fluorocarbon  liquids  in  which  the  components 
are  immersed. 

C.  STUDIES  ON  NATURAL  CONVECTION  COOLING  OF  ELECTRONIC 

DEVICES 

1.  Numerical  Studies 

Among  the  different  convection  heat  transfer  modes  mentioned  above, 
immersion  cooling  by  means  of  natural  convection  has  received  much  attention 
recently  because  of  its  high  dissipation  capabilities,  together  with  such  added 
advantages  as  no  noise  and  high  reliability  [Ref.  7],  Starting  in  early  1960  s,  numerous 
studies  of  natural  convection  in  rectangular  enclosures  have  been  carried  out.  Many  of 
the  initial  studies  were  however  motivated  by  applications  such  as  energy  transfer 
within  buildings,  solar  collection  etc..  Only  recently  has  computational  attention  been 
focused  on  simulations  of  natural  convection  in  electronic  component  cooling. 

Acharya  et.  al  [Ref.  8)  have  investigated  the  natural  convection  phenomenon 
in  an  inclined  square  box  shaped  enclosure  which  has  energy  sources  in  it.  Torok 
[Ref.  9j  used  finite  elements  technique  to  predict  component  and  board  temperatures 
and  heat  fluxes.  Kuhn  et.  al  [Ref  10]  have  studied  the  three  dimensional  natural 
convection  heat  transfer  in  a  rectangular  enclosure  with  heating  elements  on  the 
vertieail  wall  inside  the  enclosure.  Thev  found  that  when  two  elements  are 


V. 

s 

t! 


5WS 


symmetrically  placed  on  the  wall,  the  distance  between  them  has  only  a  little  effect  on 
the  mean  heat  transfer  rate  distribution,  but  has  a  significant  effect  on  the  local  Nusscit 
Numbers.  Shakerin  et.  al  [Ref.  11]  have  investigated  numerically  the  effect  of  the 
discrete  roughness  elements  protruding  from  the  wall  of  an  enclosure.  They  repeated 
the  numerical  procedure  for  both  a  single  protruding  element  and  two  elements. 

Liu  et.  al  [Ref.  7]  have  carried  out  the  finite  difference  numerical  study  of 
natural  convection  in  a  rectangular  enclosure  with  an  array  of  chips  mounted  on  a 
vertical  wall.  The  side  walls  of  the  enclosure  are  assumed  to  be  adiabatic.  The  result  of 
this  numerical  study  have  been  presented  as  temperatures  and  velocity  fields  for  the 
enclosure  widths  of  18  and  30  mm.  It  was  concluded  that: 

•  Temperature  field  in  the  enclosure  is  characterized  by  boundary  layer  regions 
surrounding  individual  chips. 

•  There  is  only  small  interference  among  the  chips,  especially  for  the  lower  rows. 

•  The  maximum  temperatures  on  the  chip  surfaces  are  located  on  the  upper 
horizontal  faces  of  the  chips.  Lesser  temperatures  have  been  found  on  the  lower 
horizontal  faces  of  the  chips  in  the  bottom  row. 

•  The  results  for  30  mm  width  are  essentially  the  same  as  those  for  IS  mm., 
indicating  only  a  minor  dependence  on  enclosure  width. 

Liu  et.  al  have  also  studied  the  local  oscillatory  surface  temperature  responses  in  the 
same  anclosure  [Ref.  12].  This  numerical  study  has  shown  that  the  maximum 
temperatures  of  chip  surfaces  have  a  tendency  to  fluctuate  within  a  range  up  to  almost 
±  3  °C  with  a  period  of  4  seconds. 

2.  Experimental  Studies 

Numerous  experimental  works  relating  to  the  enclosure  flows  have  been 
conducted  along  with  the  numerical  studies  mentioned  above.  Park  et.  al  [Ref.  13]  used 
thin  foil  heaters  to  simulate  both  flush  mounted  and  protruding  microelectronic 
components.  They  obtained  heat  transfer  coefficients  for  various  heater  heights  and 
widths.  Knock  [Ref.  2]  studied  the  effect  of  the  location  of  a  single  protruding  heater  in 
an  enclosure,  using  water  as  the  working  fluid.  Experiments  were  conducted  with  the 
heater  near  the  bottom,  at  the  center  and  near  the  top  of  a  vertical  side  wall.  Side  walls 
were  insulated  to  insure  adiabatic  boundary  conditions  on  them.  He  repeated  the 
experimental  procedure  for  three  heater  locations.  The  heater  position  was  varied  from 
close  to  bottom  of  a  side  wall  to  near  the  top  of  the  side  wall.  Mis  resuits  can  be 
quantitatively  summarized  as  a  trend  suggesting  that  as  the  heater  is  raised  within  the 
enclosure,  the  N'usselt  number  decreases. 


Filis  et.  a!  [Ref.  14]  have  investigated  the  effect  of  wall  temperature 
nonuniformity  on  natural  convection  in  an  enclosure.  They  have  presented  their  results 
stressing  that  the  nonuniform  wall  temperature  is  a  more  realistic  boundary  condition. 
Anderson  et.  al  [Ref.  15]  have  studied  the  wail  roughness  effect  on  heat  transfer 
enhancement  in  a  water  filled  cubical  enclosure. 

D.  OBJECTIVES 

The  experimental  work  described  here  was  motivated  by  the  numerical  study  by- 
Liu  et  al.  [Ref.  7]  and  is  the  basis  of  the  experimental  work  of  this  thesis.  The 
objectives  of  this  study  are  : 

•  To  design  and  build  an  enclosure  to  cool  an  array  of  nine  simulated  chips 
immersed  in  FC-75.  These  simulated  chips  are  to  be  instrumented  with  a 
number  of  thermocouples,  which  are  to  be  monitored  using  a  data  acquisition 
system. 

•  To  analyze  the  data  in  terms  of  relevant  dimensionless  heat  transfer  parameters 
such  as  Rayleigh  vs.  Nussclt  number. 

•  To  prepare  a  basis  for  similar  future  experiments,  including  various  chip 
configurations  and  enclosure  widths. 
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II.  EXPERIMENTAL  APPARATUS 

A.  GENERAL  CONSIDERATIONS 

The  apparatus  used  in  the  experiment  consists  of  an  enclosure  filled  with 
dialectic  liquid.  The  top  and  bottom  surfaces  of  the  enclosure  are  formed  by  two  heat 
exchangers  maintained  at  prescribed  temperatures  through  the  use  of  a  water 
circulation  bath.  Simulated  electronic  components  are  mounted  on  one  of  the  vertical 
sidewalls.  All  other  walls  are  exposed  to  ambient  air.  Top  view  of  the  enclosure  is  seen 
in  Figure  2.1. 

A  simulated  circuit  board  seen  in  Figure  2.2  contains  a  3  by  3  array  of 
symmetrically  placed  heating  elements.  The  heating  element  dimensions  correspond 
approximately  to  those  of  a  20-pin  Dual-in-line-Package  (DIP).  The  dimensions  of  the 
heating  element  used  in  the  present  study  were  identical  to  those  in  [Ref.  7]  and 
[Ref.  12]  to  allow  future  comparisons  between  experiments  and  computations.  The 
heating  element  was  a  rectangular  parallelopiped  made  of  aluminum  with  a  foil  typed 
heater  attached  to  the  bottom  face. 

The  present  study  dealt  with  temperature  measurements  on  the  surfaces  of 
various  heating  elements,  for  a  number  of  different  power  levels.  The  temperatures 
were  measured  by  fine  thermocouples  attaced  to  each  aluminum  block  surfaces.  The 
data  was  acquired  using  a  micro-computer  controlled  data-acquisition  system.  Details 
of  the  major  elements  of  the  experimental  assembly  are  described  next. 

B.  COMPONENTS 

1.  Heating  Element 

The  heating  element  is  6  mm  wide,  8  mm  long  and  24  mm  high.  (Figure  2.3). 
In  order  to  be  able  to  measure  the  temperatures  of  the  surfaces,  the  lateral  surfaces  are 
grooved  to  a  0.5  mm  depth  to  contain  thermocouple  beads.  The  front  surface  of  the 
chip  is  reached  by  a  hole  drilled  through  the  thickness  of  the  component.  The  beads 
and  other  uninsulated  parts  of  the  thermocouples  are  insulated  using  a  thin  layer  of 
electrical  insulating  varnish  to  insure  a  high  enough  electrical  resistance  between  the 
beads  and  aluminum  block.  An  electrical  contact  described  above  may  cause  extra 
bimetallic  junctions,  resulting  in  false  temperature  read-outs.  The  rest  of  the  groove  left 
from  thermocouple  bead  was  filled  with  a  high  thermal  conductivity  epoxy,  OMEGA 
BOND  101.  The  surfaces  of  the  blocks  were  smoothed  out  before  the  epoxy  cured. 
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2.  Thermocouples 

Because  of  the  relatively  small  size  of  the  aluminum  blocks,  a  \ery  fine.  3-nul 
copper-constantan  (Omega)  thermocouple  wire  was  chosen.  By  doing  so.  only  a  small 
amount  of  alununum  had  to  be  drilled  out  of  the  enure  block  (Approximately  "  turn3 
vs.  1152  inm3,  or  in  other  words  0.6 %  m  volume  or  weight.!.  Suice  it  is  quite  Jiiliuilt 
to  handle  these  thermocouple  wires,  special  care  had  to  be  taken  when  making  the 
thermocouples.  On  the  thermocouple  welder,  current  and  arc  time  should  be  set  to 
minimum'  and  the  argon  gas  pressure  should  not  exceed  5  psig.  The  copper  wire 
holder  of  the  welder  should  be  cleaned  of  carbon  build-up  with  a  fine  sand  paper 
frequently  to  maintain  the  high  electrical  conductivity  to  the  arc  path,  tacn  bead 


should  be  checked  for  strength  before  using,  since  after  the  cooling  oil.  the  bead 
junction  might  get  brittle  and  break  even  if  it  looks  welded. 

3.  Heaters 

Marchi  Associates  foil  type  heaters  were  used  to  provide  a  nearly  constant 
heat  flux  boundary  condition  at  the  back  cf  the  block  (Figure  2.4).  Each  of  the 
resistances  were  measured  to  be  between  10.5  to  [1.8  Q's.  The  resistor  loop  conductor 
material  is  Inconel-600.  The  back  of  the  heater  element  is  laminated  with  Kapton  with 
a  maximum  allowable  operating  temperature  of  about  250°C.  The  thickness  of  the 
heater  element  is  0.18  mm.  Each  of  the  heater  elements  has  holes  at  the  same  locations 
as  the  back  of  the  aluminum  blocks  to  pass  the  thermocouple  wire  terminating  on  the 
front  face.  In  addition,  the  heating  element  has  four  slots  on  the  sides  to  allow  passage 
of  other  thermocouples. 

The  power  leads  from  the  heater  resistor  loops  are  gold  plated  to  insure  low 
electrical  contact  resistance  on  soldering.  Wires  are  soldered  to  gold  plated  power  lead 
hook-ups.  A  length  of  5  mm  from  the  bottom  of  the  foil  heaters  had  to  be  clipped  olT 
in  order  to  accomodate  a  6  mm  distance  between  blocks.  Heaters  were  bonded  to  the 
blocks  using  OMEGA  BOND  100  adhesive  epoxy,  after  making  sure  that  the  holes  and 
blocks  were  lined  up  by  means  of  a  piece  of  wire.  They  were  then  squeezed  together  in 
a  small  hand  vice.  Teflon  sheets  were  used  on  both  faces  of  the  vice  to  protect  the 
surfaces  of  the  heater  and  aluminum  block. 

4.  Simulated  Circuit  Board 

The  vertical  sidewalls  containing  the  heating  element  arc  made  of  a  5.5  mm 
thick  plexiglass  sheet.  The  120  mm  by  144  mm  card  was  grooved  at  the  heater 
locations  to  a  depth  of  0.5  mm  to  eliminate  the  extra  thickness  that  would  have 
resulted  from  the  adhesive  used  to  bond  both  heater  to  simulated  chip  and  chip  to 
circuit  board.  The  components  with  heaters  were  bonded  to  the  card  using  rubber  to 
metal  cement.  After  placing  the  thermocouples  and  clipping  oil  the  extra  length  on  the 
foil  heaters,  the  grooves  were  filled  with  OMEGA  BOND  101  fast  set  epoxy  and 
smoothed  out. 

5.  Test  Chamber 

The  test  chamber  is  made  of  13  mm  thick  plexiglass  which  has  the  dimensions 
of  120  mm  length,  144  mm  height  and  30  mm  width.  The  spacing  at  the  back  is 
provided  to  earn.’  the  wiring  out  of  the  test  section.  All  the  wiring  was  taken  out 
through  a  Tygon  guide  tube.  Two  3  mm  aluminum  plates  were  used  as  lids  for  the 
bottom  and  the  top  of  the  enclosure. 


igure  2.6  Photograph  of  the  Test  Chamber. 


6,  Heat  Exchangers 

The  two  aluminum  heat  exchangers  measure  38  mm  by  65  mm  by  274  mm. 
The  details  of  the  heat  exchangers  are  given  in  [Ref.  2].  The  top  and  bottom  aluminum 
lids  were  grooved  to  a  2.5  mm  depth  to  provide  thermocouple  slots.  This  left  a  0.5  nun 
wall  thickness  of  aluminum  between  the  thermocouple  bead  and  the  other  surface  of 
the  aluminum  plate.  The  thermal  resistance  of  the  0.5  mm  aluminum  wall  was 
neglected.  Heat  exchangers  were  fed  with  water  as  coolant  (in  series)  by  means  of 
tygon  tubing.  Four  10-mil  thermocouples  were  used  to  read  the  top  and  bottom  heat 
exchanger  temperatures. 

7.  Assembly 

After  sliding  the  card  into  the  slot  in  the  enclosure,  both  thermocouple  and 
power  wires  were  run  through  a  small  1  2  in.  diameter  extension  pipe  behind  the 
enclosure  and  then  taken  out  through  a  1  2  in.  diameter  tygon  tube  of  25  cm.  length. 
Top  and  bottom  heat  exchangers  were  bonded  to  the  encloure  using  RTV  732 
adhesive,  sealant.  The  back  space  of  the  enclosure  was  filled  with  rubber  pieces  to 
decrease  conduction  losses  to  the  back  of  the  card.  After  attaching  the  heat  exchangers 
to  the  enclosure,  the  bottom  heat  exchanger  was  fixed  on  an  adjustable  level  aluminum 
stand  by  means  of  two  aluminum  clamps.  A  4  cm.  thick  rubber  slab  was  placed 
between  the  aluminum  stand  and  the  bottom  heat  exchanger  to  absorb  the  vibrations 
from  the  ground.  Finally,  inlet  and  outlet  tygon  tubes  from  the  constant  temperature 
bath  were  connected  to  the  heat  exchangers.  The  level  of  the  apparatus  was  adjusted 
by  means  of  the  adjustment  screws  under  the  aluminum  stand. 

C.  INSTRUMENTATION 
i .  Power  to  the  Chips 

Each  of  the  nine  heaters  was  connected  in  series  with  a  seperate  2  Ohm 
precision  resistor.  The  nine  elements  were  connected  in  parallel  to  the  power  supply. 
Power  input  to  each  heating  element  was  calculated  from  a  simultaneous  measurement 
of  the  voltage  drop  across  the  precision  resistor  and  the  overall  voltage  drop,  (see 
Appendix  A) 

The  power  panel  has  one  main  input  and  10  terminals,  the  last  one  being  a 
spare  terminal  in  this  study.  The  main  input  is  fed  from  a  DC  power  supply.  (KEPCO 
0-100  V,  0-5  A,  Mod.  JQE  100-5  type) 


2.  Data  Acquisition  System 

A  Hewlett-Packard  3054A  Automatic  Data  Acquisition  Control  System  was 
used  to  collect  the  data.  The  first  49  channels  of  the  data  acquisition  system  are  for  the 
thermocouple  readings.  Of  these,  the  first  forty-five  are  used  for  temperature 
measurements  of  the  heating  elements.  Channels  46-49  are  for  the  heat  exchanger 
temperature  measurements.  Channel  50  reads  the  main  input  voltage  to  the  power 
panel.  Channels  51-59  read  individual  heating  element  voltages. 

3.  Computer 

A  Hewlett-Packard  9826  computer  with  Hewlett-Packard  2671A  printer  was 
used  as  the  controller  for  the  data  acquisition  system.  Two  programs  written  in  HP- 
BASIC  language  have  the  data  acquisition  system  collect  the  data  and  store  in  a  data 
file  for  each  run. 


III.  EXPERIMENTAL  PROCEDURE 


A.  APPARATUS  PREPERATION 

After  assembling  the  apparatus,  the  enclosure  was  filled  through  the  wire  lead 
tube  folded  90°  upward,  using  a  small  funnel.  To  damp  disturbances  from  the 
surroundings,  the  apparatus  is  placed  on  a  seperate  table.  The  bubbles  inside  the 
enclosure  were  driven  out  by  laying  the  enclosure  horizontally  and  tilting  it  around. 
After  insuring  that  the  liquid  is  in  contact  with  the  top  heat  exchanger,  the  enclosure 
was  fixed  on  the  aluminum  stand  using  the  aluminum  clamps,  with  a  rubber  vibration 
damper  in  between.  The  apparatus  was  then  leveled  by  means  of  a  bu'obie  level  placed 
on  the  top  heat  exchanger.  Tygon  tubing  feeding  coolant  water  to  the  heat  exchangers 
were  attached  to  them  in  series.  In  other  words,  outlet  of  the  top  heat  exchanger  was 
the  inlet  of  the  bottom  one.  Because  of  the  very  high  coolant  mass  flow  rate,  the 
temperature  drop  through  the  heat  exchangers  was  negligible.  Ten  terminals  on  the 
power  panel  were  used  to  supply  power  to  the  heaters  on  the  simulated  chips.  The  first 
terminal  was  the  main  input,  the  next  nine  are  the  terminals  for  each  heater.  They  were 
connected  in  parallel  to  the  main  input  with  2CI  resistors  which  were  in  series  with 
each  heater. 

B.  DATA  ACQUISITION 

All  the  thermocouples  were  initially  scanned  to  insure  that  they  were  reading  the 
same  temperature  at  steady  state,  for  a  predetermined  sink  temperature.  The  heat 
exchangers  were  brought  to  21°C  by  means  of  the  adjustment  dial  on  the  constant 
temperature  bath.  The  voltage  on  the  power  supply  was  next  set  to  the  desired  level. 
The  voltage  setting  was  calculated  knowing  that  each  heater  and  resistor  have  a  total 
resistance  of  approximately  13STs.  Once  the  desired  power  was  determined,  voltage  on 
the  power  supply  could  be  set  to  its  approximate  value. 

Once  the  voltage  on  the  power  supply  and  the  temperature  on  the  constant 
temperature  bath  was  set  to  their  desired  values,  all  the  readings  were  taken  by  running 
the  data  acquisition  program. 

Maximum  power  per  component  was  about  2.35  Watts.  At  higher  values, 
degasing  of  the  bonding  material  resulted  in  bubble  formation  at  the  heater  blocks. 
For  each  seperate  power  setting,  sescral  runs  with  1  hour  time  period  were  taken. 


When  the  temperature  change  between  two  successive  readings  was  within  2%,  steady 
conditions  were  assumed.  This  process  takes  2-4  hours  for  each  run,  depending  on  the 
power  ievei. 

C.  DATA  ANALYSIS 

For  each  seperate  run,  all  data  were  stored  in  files  created  by  the  data  acquisition 
program.  After  the  end  of  all  the  runs,  steady  state  data  were  kept  on  the  disk.  These 
were  accessed  by  a  program  written  for  the  heat  transfer  calculations. 

All  the  raw  data  were  used  for  determining  Nusselt  and  Rayleigh  numbers  for 
each  chip  at  each  power  level. 

1.  Determination  of  Nusselt  Number 

Nusselt  Number  is  defined  as  following  : 

h  L 

Xu  =  -  (eqn  3.1) 

kf 

Where  L  is  the  characteristic  length,  taken  here  as  the  component  height.  All  the  fiuid 
properties  were  taken  variable  with  average  film  temperature.  In  order  to  evaluate  the 
heat  transfer  coefficient  h,  an  estimate  for  conduction  loss  through  the  back  of  the  test 
surface  has  to  be  made.  A  series  of  conduction  thermal  resistances,  together  with  a 
convection  boundary  condition  at  the  outside  wall  of  the  enclosure  determines  the  total 
heat  resistance.  (Appendix  A  ) 

Conduction  resistance  Rcon  for  each  solid  layer  is  determined  from  : 


L 


k  A 


(eqn  3.2) 


Where  L  is  the  length  along  the  conduction  path,  k  is  the  thermal  conductivity  of  the 
material  (  see  Appendix  A  )  and  A  is  the  area  normal  to  the  conduction  path.  Then 
the  conduction  loss  is : 


K» 


T,-Ta 


Q 

^con 


Power  input  to  each  chip  is  calculated  from 


(eqn  3.3) 


(  V.  -VJV. 

'  in  h '  h 


Qin  = 


(eqn  3.4) 


Where  V.n  is  the  input  voltage  applied  to  both  2Q  resistor  and  the  heater  which  is  in 

series  with  it.  V.  is  the  voltaee  drop  across  the  heater.  The  difference  between  V.  and 
>1  in 

V.  gives  the  voltage  crop  on  the  2fi  resistor. 

The  net  energy  added  to  the  fluid  in  the  enclosure  is  the  difference  between 
Qln  and  Qcon  which  is  at  the  same  time  from  Newton's  law  of  cooling  : 

Qconv=  kA(Tavg“Tc) 

-  <3;n-Q«=„ 


where. 


A  =  Total  surface  area  of  the  chip. 

T  *  Area  weighted  average  surface  temperature. 
T  =  average  heat  exchanger  temperature. 


I  Vi 


T 

avg 


(eqn  3.5) 
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Tj  =  temperature  of  i*  surface  of  the  chip, 
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Aj  =  area  of  the  1th  surface  of  the  chip.  The  summation 


index  i.  varies  from  1  to  5. 


The  above  procedure  is  repeated  for  each  of  the  nine  chips  at  each  power 
level.  Then,  from  Equation  3.1,  Nusselt  number  is  calculated. 

2.  Determination  of  Rayleigh  number 
The  definition  of  Ravleiah  number  is 


Ra  =  (Gri(Pr) 


Where 


Gr  = 


(eqn  3.6) 


Pr  =  v  a 
a  =  kf(pcp) 

All  the  properties  in  the  above  equations  are  calculated  at  average  film 


temperature 


T  +  T 

avg  c 


(eqn  3.7') 


Functional  relationships  for  the  properties  at  the  film  temperatures  are  given 
in  Appendix  A,  together  with  the  sample  calculations.  Uncertainty  calculations  are 
given  in  Appendix  B. 

3.  Flux  Based  Rayleigh  Number 

Another  approach  is  to  calculate  a  modified  Rayleigh  Number  based  on  the 
heat  flux.  This  Rayleigh  number  is  defined  as: 


Ra  =  (Gr  )(Pr) 

me  me  '  ’ 

Where  "c"  stands  for  the  center  of  the  heated  block. 
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(eqn  3.S) 


where  q"  =  Q  A  is  the  heat  flux  at  the  surface  of  the  chip  and  ; 


Nuc  = 


(eqn  3.9) 


h  =  q'V  (Tavg-Tf) 

The  above  definition  of  Ramc  and  Nuc  have  been  used  to  compare  the  results 
with  those  of  Fujii  and  Fujii  [Ref.  23].  Fujii  and  Fujii  correlation  of  boundary  layer 
solution  for  a  vertical  surface  with  constant  heat  flux  is  given  by  : 


(eqn  3.10) 


where 


Nuc  =  fIPr)(GrmcPr)° 


fljPr)  =  (Pr/(4  +  9VPr+  10Pr))°'2 
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IV.  RESULTS 


A.  TEMPERATURE  DATA 

As  indicated  in  Chapter  III,  the  raw  temperature  data  were  kept  in  order  to  be 
processed  later  on  by  the  heat  transfer  calculations  program.  Temperature  outputs 
belonging  to  6  power  levels  are  as  following  (see  figure  2.2  for  identification  of  chip 
numbers). 


TABLE  1 

TEMPERATURE  DATA  FOR  INPUT  POWER  Q,N  =  0.34  W 


DATA  IS 

STORED 

IN  DATA 

FILE 

:  DATA7 

CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

POWER(WATTS) 

CHIP 

N01 : 

29.4 

29.7 

29.7 

28.9 

29.8 

.36 

CHIP 

N02: 

31 .4 

31 .6 

30.9 

29.7 

31 . 1 

.35 

CHIP 

N03: 

31 .9 

32.0 

31 .9 

31.9 

32.0 

.33 

CHIP 

N04: 

30.0 

28.6 

29.2 

28.4 

29.4 

.35 

CHIP 

NOS: 

31  .8 

31  .6 

31.1 

31.7 

31.0 

.35 

CHIP 

N06  : 

32.3 

32. S 

31.4 

31.6 

32.5 

.35 

CHIP 

N07 : 

30.1 

32.0 

29.4 

29.4 

32.0 

.35 

CHIP 

N08: 

31 .2 

31 . 1 

30.  S 

30.8 

31  .3 

.32 

CHIP 

N09: 

32. S 

31  .S 

31.3 

31  .3 

32.0 

.36 

HEAT  EXCHANGER  TEMPERATURES:  RIGHT  LEFT 

BOT :  20.7  20.6 

TOP:  20.8  20.7 
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TABLE  2 
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TEMPERATURE  DATA  FOR  INPUT  POWER  Q,N  =  0.68  W 


DATA  IS 

STORED 

IN  DATA 

FILE 

:  DATA  12 

CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

POWER  <  WATTS) 

CHIP 

NOt : 

34.2 

34.8 

34.8 

33 . 5 

35.0 

.68 

CHIP 

N02: 

37.9 

37.5 

37.0 

34.9 

37.3 

.69 

CHIP 

N03: 

38.9 

39.3 

39.1 

39.1 

39.2 

.  66 

CHIP 

N04: 

35.4 

33.0 

34.0 

34.3 

34.3 

.68 

CHIP 

NOS: 

38.4 

38.0 

37.3 

38.2 

37.0 

.  68 

CHIP 

N06: 

39.3 

39.8 

37.4 

38.0 

39.8 

.69 

CHIP 

N07 : 

35.6 

35.1 

34.2 

34.2 

35.4 

.68 

CHIP 

N08: 

37.5 

37.3 

36.4 

36.8 

37.8 

.64 

CHIP 

NO  9: 

39.8 

38.0 

37.2 

37.2 

39.2 

.68 

HEAT  EXCHANGER  TEMPERATURES:  RIGHT  LEFT 

BOT:  20.8  20.5 

TOP:  20.9  20. S 
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TABLE  3 

TEMPERATURE  DATA  FOR  INPUT  POWER  Q,N  =  1.07  W 


DATA  IS 

STORED 

IN  DATA 

FILE 

:  DATA16 

CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

POWER  (Uf 

CHIP 

N01 : 

40.9 

41  .S 

41  .5 

39.7 

42.2 

1 .08 

CHIP 

NO  2: 

46 . 5 

4S.6 

4S.2 

42.3 

46.0 

1.09 

CHIP 

N03: 

48.1 

48.9 

48. S 

48. S 

48.7 

1  .OS 

CHIP 

N04: 

43.0 

38.8 

40.9 

41  .3 

41  .3 

1  .08 

CHIP 

NOS: 

47. S 

46.4 

4S.7 

47.1 

4S.2 

1  .08 

CHIP 

N06: 

48.6 

49.2 

4S.S 

46.2 

49.3 

1.09 

CHIP 

NO  7  : 

43.2 

41  .9 

40.9 

40.9 

42.6 

1  .08 

CHIP 

N08: 

46.2 

4S.S 

44. S 

4S .  1 

46.6 

1.02 

CHIP 

NOS: 

4S.4 

47.1 

4S.6 

4S.6 

48.7 

1  .09 

HEAT  EXCHANGER  TEMPERATURES:  RIGHT  LEFT 

BOT:  21.1  20.6 

TOP:  21.2  20.7 
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TABLE  5 

TEMPERATURE  DATA  FOR  INPUT  POWER  Q..  =  1.95  W 


DATA  IS 

STORED 

IN  DATA 

FILE 

:  DATA27 

CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

P0WER<  Wf 

CHIP 

NOT : 

SS.1 

S5.7 

56.2 

52.8 

57.6 

1  .95 

CHIP 

NO  2: 

63.  S 

62.4 

62.3 

59.2 

62.5 

1.97 

CHIP 

N03: 

66 . 3 

77.2 

71  .8 

72.5 

74.8 

1  .94 

CHIP 

NO  4: 

S9 . 1 

51.3 

55.5 

56 . 3 

56.0 

1  .95 

CHIP 

NOS: 

6S.S 

62.2 

62.8 

64.8 

61.7 

1  .96 

CHIP 

NO  6 : 

66.0 

67. S 

61 .9 

62.8 

67.3 

1  .97 

CHIP 

N07: 

S9.S 

56. 2 

55.7 

56.2 

58.1 

1  .95 

CHIP 

N08: 

63.8 

62.4 

60.2 

61.9 

64.3 

1.88 

CHIP 

N09: 

68.2 

64.7 

63.1 

63.7 

75.5 

1  .97 

HEAT  EXCHANGER  TEMPERATURES:  RIGHT  LEFT 

BOT:  22.2  20.8 

TOP:  22.4  20.9 
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TABLE  6 

TEMPERATURE  DATA  FOR  INPUT  POWER  Q,s  =  2.35  W 


DATA  IS 

STORED 

IN  DATA 

FILE 

:  DATA32 

CENTER 

TOP 

RIGHT 

LEFT 

BOTTOM 

P0WER( WATTS ) 

CHIP 

N01  : 

64.  9 

64.9 

bS.S 

61  .8 

67.4 

2.3S 

t 

CHIP 

N02  : 

73.2 

71 .8 

71  .6 

68 .2 

72.0 

2.38 

CHIP 

N03: 

7S.6 

83.3 

79. S 

80.2 

81 .7 

2.3S 

CHIP 

N04: 

68.6 

60.4 

64.8 

6*4.9 

63.2 

2.36 

CHIP 

NOS: 

7S.0 

71  .5 

72.1 

74.9 

71  .7 

2.37 

CHIP 

N06 : 

7S.6 

76.7 

70.7 

72.4 

76.6 

2.38 

CHIP 

N07 : 

70.3 

6S.S 

63 . 3 

66.1 

o  8 . 0 

-  ->C 

•_  •  w  -* 

CHIP 

NOB: 

73.9 

72.6 

70.3 

71.8 

74.  S 

2.26 

CHIP 

N09: 

77.9 

72. S 

72.9 

72.6 

'■  z 

Ow  •  -/ 

2.37 

HEAT  EXCHANGER  TEMPERATURES:  RIGHT  LEFT 

BOT:  23.8  21.0 

TOP:  24.1  2:. 2 
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Following  are  the  graphical  representation  of  the  temperature  data  given  above. 
The  plots  are  prepared  to  include  the  center,  maximum  and  minimum  temperatures  for 
each  chip  at  6  power  levels. 
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B.  HEAT  TRANSFER  RESULTS 

The  temperature  data  given  in  the  previous  section  have  been  processed  by  the 
heat  transfer  calculations  program  given  in  Appendix  C.  Following  are  the  tabulation 
of  these  heat  transfer  calculations  : 


If 


% 

I. 

9 

l 

| 

i 

1 

* 


I 


| 

I 


I 


QNET(W) 
.  33S 


.324 

.292 


.320 


.280 


QNET(W) 

.66  2 

.600 


.574 

.636 

.570 

.564 


ft 

1 

£ 


TABLE  ’ 

HEAT  TRANSFER  CALC.L  LATIONS  I  OR.  0IN  =  '>.34  W 


TAVG-TS 

8.72 

10.17 

11.26 


1  0 . 34 

11.25 
3.38 

10.26 
11.08 


THE  RAW  EMF  DATA  ARE  FROM  FILE:  DATA7 


~>c  op 

w  -J  •  JO 


19.49 
1  6 . 32 


24.37 


17.81 

16.73 


c  c 


6.75 


XUNC.IN  NU  RA+1.E-7  7.UNC 


26.78 
30.73 
33.33 
27.30 
31  .  76 
32.97 


28.1  8 


5.59 
6 . 62 


C  -  o 
-•  •  o 


7.10 


7 . 4U 


6 . 06 


■|  cl 

i  .  _  o 


IN  RA 

89 


TABLE  S 


HEAT  TRANSFER  CALCULATIONS  TOR  Q,s  -  O.oS  \V 


TAVG-TS 

13.58 

16.16 

18.38 
13.83 
17.26 

17.38 
14. 17 

16.38 
17.63 


THE  RAW  EMF  DATA  ARE  FROM  FILE:  DATA12 


NU 

32.33 

24.67 

20.22 

30.67 


21.35 

29.75 

23.12 


21  .26 


7.UNC.IN  NU  RA*1.E-7  7.UNC.IN  RA 


22.10 

25.09 

27.18 

23.01 

25.82 

26.73 

23.16 

26.1 1 

27  .'44 


9.15 

11.16 

12.98 

9.34 

12.06 

12.57 

9.60 

11.34 

12.36 


2.85 

2.91 


2.85 

2.95 

2.97 


2.92 

2.96 
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TABLE  9 

HEAT  TRANSFER  CALCULATIONS  I  OR  Q,s  =  I W 


THE  RAW  EMF  DATA  ARE  FROM  FILE:  DATA16 


QNET(W) 

1.054 

TAVG-TS 

19.98 

NU 

35.14 

7.UNC.IN  NU 
20.93 

RA* 1 . E- 7 
14.40 

7.UNC.  IN 

3.1)7 

.960 

24.  t  I 

26.60 

23.65 

18.07 

3.25 

.898 

27.54 

21  .83 

25.50 

21  .32 

3.41 

1.016 

20.64 

1°  70 

•  1  J 

21  .36 

14.96 

.  1 0 

.938 

25.78 

"*/*  A  “7 

24 . 45 

19.63 

“»  *)  ’* 

.921 

26 . 42 

23.3! 

25.13 

20.24 

3.35 

1.010 

20.98 

32 . 0  7 

22.01 

15.27 

3.1 1 

.908 

24.5  b 

24. 70 

24.75 

1  6 . 49 

~  A  7 

.906 

26 . 32 

->  :■  o  '■ 

25.74 

20.14 

:  ->c 

TABLE  10 

HEAT  TRANSFER  CALCULATIONS  FOR  Q,N  =  Leo  \V 

THE  RAW  EMF  DATA  ARE  FROM  FILE:  DATA22 


QNET(W) 

1  .61 1 

TAVG-TS 

28.99 

NU 

37.23 

XUNC.IN  NU 
20.02 

RA-*  1  .  E  -  7 

22.80 

7.UNC.  IN 
3.49 

1.486 

35.11 

28.47 

22.28 

29.20 

3.79 

1  .416 

40.61 

23.54 

23.65 

35 . 48 

4. 08 

1  .547 

30.37 

34 .  1 6 

21.01 

24.18 

3.55 

1.458 

36.39 

26.96 

22.84 

30.62 

3.86 

1.428 

37.97 

25. 34 

23 . 62 

32.41 

3.34 

1  .540 

30.72 

*»  r,  - 
Jj  .  Ow 

21.13 

24 . 54 

3.57 

1.420 

55.07 

^  7  .  ^  ’’ 

■*  *  * 

29.  16 

3.79 

1  .406 

38 . 1  i 

24.86 

24.19 

32.57 

3.95 

« 

fa 
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TABLE  II 

HEAT  TRANSFER  CALCULATIONS  TOR  Q,s  =  1.95  \V 


THE  RAW  EMF  DATA  ARE  FROM  FILE:  DATA27 


ONET(W) 

1 .897 

TAVG-TS 

33.49 

NU 

38.08 

7.UNC.IN  NU 
19.87 

RA-*  1  .  E-  7 
27.73 

7.UNC.  IN 
3.77 

1.759 

40.40 

29.39 

21.39 

35.57 

4.13 

1  .686 

49.27 

23.23 

}  o  ^  t» 

46.84 

4.62 

1.826 

35.06 

35 . 05 

20.82 

25.43 

3.85 

1.714 

42.53 

27.26 

22.7  7 

38 . 1 5 

4.24 

1.713 

42.35 

27 . 04 

22.86 

38.54 

4.26 

1.817 

35.75 

84.22 

20.94 

30.21 

3.38 

1.681 

40.31 

*1  "7  •-*  A 

w  .  •  O  L- 

22 . 30 

8  6 . 06 

4.15 

1  .673 

44.56 

A  c  •> 

4- 

22.7  0 

40.69 

4.36 

TABLE  12 

HEA  L  TRANSFER  CALCULATIONS  TOR  Q|N  =  2.35  \V 


THE  RAW  EMF  DATA  ARE  FROM  FILE:  DATA32 


•M 


ONET(W) 

2.257 

TAVG-TS 

42.01 

NU 

36.37 

7.UNC.IN  Nil 
20.41 

RA* 1 . E- 7 
38.11 

y.  unc.  in 

4.32 

2.125 

48.93 

29 . 52 

22.17 

47.09 

4.70 

2.062 

56.40 

24.97 

23.04 

57.76 

5.14 

2.192 

43.29 

34.31 

21.16 

39.71 

4.39 

2.082 

51  .23 

27.67 

50.27 

4.84 

2.080 

51  .27 

27.63 

22.86 

50.32 

4 . 84 

2.159 

44.96 

32.56 

21  .58 

41.84 

4 . 48 

2.020 

49 . 95 

2  7.51 

23.24 

48 . 44 

<u.76 

2.037 

53.03 

26.18 

23  v  59 

62 .81 

4.34 

46 
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Following  are  the  graphical  representation  of  heat  transfer  calculations  to  include 
Nusselt  number  versus  Rayleigh  number.  The  plots  have  been  classified  as  BOTTOM 
row  of  the  chips,  MIDDLE  row  of  the  chips  and  TOP  row  of  the  chips. 
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V.  DISCUSSIONS  AND  RECOMMENDATIONS 


A.  DISCUSSIONS  OF  THE  RESULTS  J 

Figures  including  the  temperature  distribution  for  different  power  levels  show  « 

explicitly  that  the  temperatures  of  the  simulated  chips  have  a  tendency  to  increase  with  j 

the  height  inside  the  enclosure.  This  suggests  that  the  chips  in  the  bottom  row  have  the  ; 

lowest  temperatures  while  the  top  ones  have  the  highest.  This  conclusion  is  also  j 

obvious  from  the  figures  that  include  Nusselt  vs.  Rayleigh  number.  Highest  heat  \ 

transfer  rates  or  in  other  words  highest  Nusselt  numbers  are  seen  in  the  lowest  row  of 
the  chips.  Another  trend  is  that  the  increase  in  Nusselt  number  (heat  transfer  rates) 
gets  smaller  with  increasing  Rayleigh  numbers.  It  is  supposed  that,  this  results  from  the 
significant  increase  in  the  bulk  temperature  of  the  fluid.  Comparison  of  the  results  with 
Fujii  and  Fujii  correlation  (Ref.  23J.  shows  the  weak  effect  of  protrusion  on  tiie  heat 
transfer.  The  bottom  row  of  the  chips  have  a  little  higher  heat  transfer  rates  than  those 
the  correlation  gives.  Especially  the  top  row  of  the  chips  have  significantly  lower  heat  1 

transfer  coefficients.  This  is  because  the  chips  in  the  bottom  row'  have  free  boundary  j 

layer  plumes  whereas  top  row’  chips  have  boundary  layer  plumes  confined  by  those  of  j 

lower  row's.  I 

t 

) 

B.  RECOMMENDATIONS 

High  uncertainty  in  the  Nusselt  number  calculations  results  mainly  from  the 
necessity  to  make  a  rough  estimation  for  conduction  losses  to  the  back  of  the  j 

enclosure.  This  requires  that  the  enclosure  w’alls  be  very  well  insulated  to  insure  a  very  ! 

small  amount  of  conduction  loss  so  that  the  uncertainty  is  low.  j 

A  flow  visualization  technique  for  FC-75  should  be  established.  Known  ! 

techniques  for  water  such  as  die  injection,  electro-chemical  (because  of  the  requirement 
for  maintaining  the  dielectricity)  or  floating  particles  (because  of  the  requirement  for 
neutral  density)  techniques  are  not  suitable  for  FC-75.  Repetition  of  the  experiment  ; 

with  different  enclosure  widths  and  component  configurations  will  be  worthwhile  so  far  ! 

as  the  collection  of  data  sets  for  general  similar  problems  is  concerned. 


APPENDIX  A 

SAMPLE  CALCULATIONS 

Through  this  appendix,  sample  calculations  will  be  based  on  chip  number  2  with 
power  input  0.35  W. 


1.  DETERMINATION  OF  INPUT  POWER 


Figure  A.l  Electrical  Network  of  Power  Input. 


From  Equation  5.4 


o.  «  :  =  0.351  w 


«? 


I 

$ 


p: 

: 


2.  ESTIMATION  OF  CONDUCTION  LOSS 

According  to  Equation  3.3.  we  need  to  know  T  and  R  .  T  can  he  caicuL.ced 
using  the  relationship  between  Qin  and  diflcrence  between  iback  surface  temperature 
of  the  chip)  and  T  (center  temperature  of  the  chip)  obtained  from  a  previous  single 
chip  experiment: 

Q.  =  <U>S3(T  “  I  ) 

in  '  s  n?' 

Since  there  are  4  chips,  we  might  want  to  use  an  average  T  .  f  rom  Table  1. 

Ts  =  35.4°C  lor  chip  no  2  and 
average  T.  =  35.0°C. 


VW — * 

Rc 

-N'C — • 
R_ 


Figure  A. 2  Thermal  Resistance  Network  tor  Conduction  Loss. 


In  order  to  estimate  the  total  thermal  resistance  R  .  we  need  to  know  the 

con 


thermal  conductivities  of  each  material,  as  given  in  Table  13. 

Ra  was  evaluated  from: 

where  ha  =  3.9S5  \V  nr-°C  as  an  average  value,  calculated  from: 


h,  =  l.42( AT  II>0  -5  from  (Ref.  22|. 
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TABLE  13 

THERMAL  CONDUCTIVITIES  OF  THE  MATERIALS 


Material 

k  (YV  m-°C) 

Rubber 

0.03S9 

Plexiglass 

0.1421 

FC-75 

0.0640 

TABLE  14 

THERMAL  RESISTANCES  TO  CONDUCTION 


Resistance 

Material 

Value  (°C  W) 

Rc 

Plexiglass 

54.979 

Rr 

Rubber 

22.321 

R. 

FC-75 

54.253 

RP 

Plexiglass 

4.073 

R 

Air 

11.690 

Rc  was  evaluated  using  a  projected  area  of  16  *  32  mm2  instead  of  8  x  24  nun2.  Rr  and 
R[  were  evaluated  assuming  they  take  up  80%  and  20%  of  the  total  heat  flow  area, 
respectively.  Since  all  the  Rc  resistances  are  parallel  to  each  other  and,  Rf  and  R,  are 
parallel  to  each  other  such  that ; 

Reqi  =  15  814  °C  w  and, 

Req2  =  6.109 ’C/W 
were  used.  Then  the  total  resistance: 

Rrn„  =  37.686  °C.  W 
con 

•  from  Equation  3.3  : 

Qcon  =  (35.0-  20.0)  37.686  =  0.398  W. 

Since,  at  the  same  time, 

Qcon  =  (Tl-VRt 


56 


where  Tj  is  the  temperature  of  the  back  of  the  circuit  board  and, 

Rt  =  Rcon  “  Req2 
=  31.577  °C  W. 

T(  =  32.57°C.  Then,  for  chip  2, 

Q  =  (35.4-32.57)  54.979  =  0.052  W. 

xcon 

3.  CALCULATION  OF  NUSSELT  AND  RAYLEIGH  NUMBER 

From  the  definitions  given  in  Chapter  III, 

T  =  20.7°C 

C 

Tavg  =  30.87°C  and, 

Qconv  =  0.351  -0.052  =  0.299  W,  the  net  power. 
Then  from  Equation  3.3, 

h  =  0.299  (5.76  x  10“ 4)(30.87 - 20.7) 

=  51.042  W,m2-°C 

From  Equation  3.7, 

Tf  =  (30.87  +  20.70)/2  =  25.79°C 

Using  the  functional  relationships  given  in  (Ref.  7]  for  properties  of  FC-75, 

kf  =  0. 1(0.65  -  7.8947  x  10  “  4  Tf) 

=  0.0629  W,  m-°C 
p  =  1000(  1.825 -0.00246Tr) 

=  1761.6  kg  m3 

Cp  =  4180(0.241111  +  3.7037  x  10“ 4Tf) 

=  1047.8  J/kg-°C 

v  =  10“6(a0  +  ajTf  +a2Tf2  +  a3Tr3  +a4Tf4) 

=  0.8526  x  10“  6  m2  s 

where 

a0  =  1.4074 
at  =  -2.964  x  [0“2 
a2  =  3.8018  x  10  “  4 
a3  =  -2.7308  x  10” 6 
a4  =  8.1679  x  10“ 9 

P  =  0.00246  (1.825 -0.00246Tf) 

=  0.0014  (°C)“l 


Then,  from  Equation  3.1  : 


<1. 


From  Equation  3.6, 


Nu  =  (51.042K0.024)  0.0629  =  19.48 
a  =  (0.0629)  (1761. 6X1047.S)  *  3.4077  x  io~sm:  s 
Pr  =  (0.8526  x  I0~b»  (3.4077  x  io_l?)  =  25.02 

Gr  -  (9.Sh(0.0014)(0.024r,( 30.87 -  20.'0|  (,S526  x  10_6); 
=  2.651  x  106 

Ra  =  (2.651  x  106)<  25.02)  =  6.63  x  (o' 


Following  is  to  show  how  Ramc  and  Nu.  are  calculated; 
From  Table  7  , 

Qnet  =  0.299  W. 


q"  =  0.299  5.76  x  10  ~  4  =  519.1  W  m2, 


•  from  Equation  3.8, 


From  Equation  3.9, 


in  which 


From  Equation  3.10, 


Gr  =  3.70  x  io6  and. 

me 

Ra  =  9.26  x  107 

me 


Nuc  =  19.4 


h  *  102.2  W/m2-°C 


Nu  =  (0.6088K39.20I)  «  23.9 
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APPENDIX  B 

UNCERTAINTY  CALCULATIONS 


Uncertainties  in  the  experiment  were  evaluated  using  the  root  mean  square 
method  of  Kline  and  McClintock  [Ref.  21 J. 

Following  is  a  summary  of  the  uncertainties  used  in  the  experiments  : 

TABLE  15 

UNCERTAINTIES  OF  THE  VARIABUES  USED  IN  THE  EXPERIMENT 


Variable 


Uncertainty 


10“  5  m 


0.2°C 


Reference 


Resolution  of  measurement  device 


Thermocouple  reading 


[Ref  2| 
[Ref  2] 


0,06 SI 


3%  variation 


[Ref.  7] 


Uncertainty  in  thermal  resistances  is  given  by  : 


where 


SR/R  =  ((SL/L)2  +  (5k  k)2  +(5A,  A)2)0'5 


6A,'A  =  ((6Lj ,'Lj)2  +«5L2  L2)2)0 


For  Rp  , 


I 

5* 

y 

i 

! 
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h 

v 
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Lj  =  0.120  m. 

L,  =  0. 144  m. 

A  =  1.723  x  10"  2  nr 


6A  A  =  0.0001.  and  ; 


6R  R  =  0.07 


6Rp  =  (4.073X0.07)  =  0.2851  °C,W. 


ft 

5* 

» 

* 

f: 


If  similar  method  is  repeated  for  all  the  thermal  resistances  ; 


8R.„„  -  «8Rra,)2  +(5R.o!>2  +<8R„)2  +(8R.)2)0  ! 


Uncertainty  in  Q.  is  dominated  by  the  uncertainty  in  the  2Q  resistance  since 
voltage  readings  done  by  the  acquisition  system  have  a  very  low  uncertainty. 

■  8<3,„  *  0  W 

For  the  chip  chosen  for  calculations  ; 

5Qin  *  (0.03)(0.351)  -  0.01053  W 

Uncertainties  for  other  variables  are  calculated  as  following  : 


5Q.o„Q=o„=  ((SCt.-t^ct.-t,))2  +  («RC,RC)2)" 


where 


6rrs-T1)  =  ((5Ts)2  +  (5Tl)2)0 


Uncertainty  in  net  power  ; 


«Q,„„v  -  «6QJ2  +(5Q,„»)2)°'s 


\ 

| 


Uncertainty  in  h 


6h  h  -  ((60  Q  )2+  (6(T  -TV(T  -T))2)0-5 

'  "  vconv  ^conv'  '  '  avg  c;  '  avg  cn  ' 


Uncertainty  in  Xusselt  number  ; 


SNu'Nu  =  ((6h  h)2+-  (6k  k)2  +  (5L  L)~)0'5 

=  ((0.3047)2  +■  (0.040)2  +  ( 10  ~  5  0.024)2)0  5  =  0.3073 

Uncertainly  in  Rayleigh  number  is  calculated  the  same  method  : 

8Ra,Ra  =  ((8Gr.Gr)2  +  (6Pr  Pr)2)0  5 

where 

SPr  Pr  =  ((8v  v)2+  (5a  a)2)0  5 
SRa  Ra  =  0.0284 


appendix  c 

PROGRAM  LISTINGS 


1.  DATA  ACQUISITION  PROGRAM 


!0  ! F  ILE  NAME: THESIS' 

20  (EDITED  3y:  TURGAT  pamuk .LTjG . TURKISH  navy 
30  ! 

40  1 

50  'PROGRAM  FOR  FATHERING  AND  REDUCING  DATA 

60  '.FOR  49  THERMOCOUPLES  AND  VOLTAGES  FOR  ®0WER 
70  (SUPPLY. THERMOCOUPLES  '  TO  45  ARE  ON  THE 

80  (BLOCKS. 46  TO  49  ARE  ON  THE  HEAT  EXCHANGERS. 
90  ! 

100  ' 

!  10  COM  /Co/  D< 7) 

120  DIM  Em r i 59 ' . Roue r ' 9 i . y  (  49  i 
130  • 

140  'CORRELATION  f aCCR'E  "C  CONVERT  EMr  TC  DEGR- 
i 50 

160  DATA  o7i 003609 ' , 2572' . 9 . -  76 7  345 . 3 . 73025696 . 

170  DATA  -9247486539.6. 39E*- . -2. 66E13.j.94£14 
i3Q  • 

: 90  READ  D i • > 

200  Ro-2.000 

JtO  PAINTER  IS  to- 

220  BEE* 

230  INPUT  "ENTER  THE,  INP'j*  **0DE  :  G-STS .  ’  *F  ILE  "  .  In 
240  IF  !«• '  Then 
2S0  BEE3 

260  INPUT  "ENTER  name  Sr  .Did-  :  i»* 

270  PRINT  USING  *2GX.** THESE  RESULTS  ARE  FROM  DAT, 
230  BEE3 

->3q  else 

5do  BEEP 

Sio  INPUT  "ENTER  name  qp  nEh  DATA  "HE* .tW ,  1*$ 
320  PR  IN  T  USING  "20R. "DATA  IS  STORED  IN  DATA  EH; 
330  BEE3 
340  END  IF 
'■SO  1 

360  IF  In*'  THEN  ASSIGN  »F i i*  TO  Oldfiiei 

•79  IF  In-0  THEN 

330  CREATE  BDAT  NewtileS.5 

390  ASSIGN  AFiie  TO  New  i  ieS 

400  END  IF 

410  ' 

420  'NOW  THE  ACQUISITION  SYSTEM  RilLi.  READ  DATA 
430  ' 

440  BEEP 

450  BEE3 

460  IF  In-0  then 

470  OUTPUT  709: ‘AR  AFOO  ALS8" 

480  0UTPUT  722:  "F!  R1  Tl  2'  F'_;- 
490  r0R  1-0  TO  S3 
S.iO  OUTPUT  709 :  *  AS  " 

510  ENTER  722:£n»*>I) 

520  BEE3 

S  30  NEXT  I 

sao  output  «f i  i«:Em»'< 4 1 
650  ElSE 

560  ENTER  *F i ie:Em» ' » ' 


(i2 


Copy  ovailohl*  to  DTIC  docs  not 
permit  fully  legible  reproduction 


570  END  Ir 

560  OUTPUT  709:  "TD* 

530  'CONVERSION  PROM  HVOLT  TO  CELCIUS 
500  EOR  1*0  TO  43 
61  0  Sum-0 . 


620 

FOR  J-0  TO  7 

530 

Sum*5um*D ( J ) 

•Em) (  I )  ‘  J 

640 

NEXT  J 

650 

Tip  *Sun 

660 

NEXT  I 

570 

i 

680 

! POWER  INPUT 

CALCULATION 

690 

J-1 

700 

Vo i t  *£mt ( 49 ) 

710 

i 

720 

FOR  1*50  TO 

58 

730 

Pow  t  J ) *£mt (  I 

) * • Vo  1 1 -Em) (1)1 /Re 

'40 

i.  » 1 

750 

NEXT  I 

7S0 
?  y 

i 

| 

7130 

BEE0 

*3,1 

BEEC 

900 

•prints  all 

THE  COMPONENT  TEMP£RATUTE5  AND  °0WERS 

61  0 

PRIN* 

320 

PRINT 

330 

PRINT  USING 

"20*.  *“  center 

TOP  RIGHT  LEFT  BOtTQM 

840 

PRINT 

550 

PRINT 

560 

PRINT  USING 

“  '2:;. '"CHIP  N0t: 

"  .5' 3D.2.2X i .4X.22.0D.  “: 

*  1  **  1 

.  C0'J  •  i  1 

;  :;i 

_ PRINT  USING 

"T2X.**CHI6  N02: 

".5; 32.0.2X) .4X.2D.DD./*: 

T '  i?  / 

. Pow i 2  ■ 

930 

°R;NT  USING 

"  1 2X , * "CHIP  N03: 

".512D.D.2X:  .4X.2D.DD.,'*: 

!  7  1 

T (  1 4  j  , row '  3  > 

330' 

PRINT  USING 

* T2X.  "CHIP  N04: 

**.5>3D.D.2X! .4X.2D.DD./*: 

H  , 

T ( 1 9 i . r  ow ' 4 i 

900 

PRINT  USING 

"TJX.  "CHIP  NOS: 

" .  5  >  30 . 2 . 2X : .  4>: .  ;o ,  00 .  X * : 

2‘j '  * 

*  > 24  j . F  ow i 5 ■ 

?!  0 

PRINT  USING 

*  1 2X .  "CHIP  N06 : 

**.5(3D.0.2X; .4X.20.DD./*: 

**o  , 

«_  j  •  • 

T  <  29 i . °ow ' 6 ) 

PRINT  USING 

“TJX. "CHI?  N07: 

"  .5 • 3D.D .2X i . 4X .22 . 00. /* : 

2  Vi . 

T  (  34  i  .  Fowl  7 

930 

PRINT  USING 

*12X,  "CHIP  N08: 

". 5(30. D. 2Xi ,4X. 20.00./*: 

38,  . 

T ( 39 ) . Pow ( 3 ) 

3«0 

PRINT  USING 

* !  2X . "CHIP  N09: 

**  .5<  3D.D.2X » .4X.2D.DD. /* : 

43  :  . 

T  1  44  1  .Fowl  9  > 

56  0 

PRINT 

960 

PRINT 

370 

i 

980 

PRINT  USING 

"5X."HEAT  EXCHANGER  TEMPERATURES:  RIGHT 

990 

PRINT  USING 

■29X.  "B0T:"  .2X 

. 2 1 30 . 2 . 2X 1  * : T i 45 • . T ( 46 , 

1000 

PRINT  USING 

*7gx  "TQP*"  *x 

,2<  30. 0.2X ! “ : T( 47 ) , T<  48  < 

1010 

BEEP 

'020 

i 

10  30 

i 

1040 

PRINTER  IS  1 

’050 

ASS IGN  *F  i  1*, 

TO  ► 

1060 

END 

:  ( o . .  t  ( t ' .  *  ■ : ; . :  < :  • . 

T  t  £  ! .  T 1 5  ,  .  r  (  7  ;  , :  i  8  ) 
T( '0 » .  T( 1 :  ! ,  r(  T 2  •  .  Ti 
T(  '5  1  ,TH6;  .  T:  !  7  '  .  T( 
r<23  :  .  2:  )  .  T.  22  > .  T' 

T ( 25  < . T ( 26  . Ti27 i . Ti 
T; 30) ,T(3i : .T(22 ' . T( 
T ( 35 : . T  <  38 ; .1(37). T( 
Ti 40 i . Ti 41 i . T(42 ) . Ti 


Copt  available  to  DTIC  does  not 
permit  fully  legible  reproduction 


2.  HEAT  TRANSFER  CALCULATIONS  PROGRAM 


!  0 
20 
30 
40 
50 
SO 
70 
80 
90 

.  Co(9) 
100 
110 
120 
130 
140 
150 
160 
170 
130 
130 
200 
210 
220 
220 
240 
250 
260 
270 
280 
290 
200 
310 
320 
330 
340 
350 
360 
370 
33  0 
390 
4  00 
410 
420 
430 

4  «J  0 

450 
480 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 


F[lE  NttME  : T 51 52 

EDI  TED  BY  :  TuRGAY  PAMUL .  LUC . 


URLI5r  NHVY 


PROGRAM  TO  ANALYZE  THE  RAW  DATA  IN  THE  DATA 
2 ILES  OBTAINED  FROM  TEST  RUNS 


CON  /Co/  D( 7 ) 

DIN  Emt i 59 ) .Power ( 9 J . T( 49 ) , 
.  N(  9 ) . Nu ( 9 ) 


Tavgi  9  i  .Ts(9>.Ttilmi9>  .Qnet  I  9)  ,H<9> .K( 9)  .RhoO  ) 


DIN  Ray (9) ,Delt(9).Alfa(9).Pr(9>.Gr(9) 

DIN  Be ta ( 9  1  . Doou (9)  . D ts ( 9 1  . Da  loss ( 9  t 
1  CORRELATION  FAC7GRS  TO  CONVERT  ENF  TO  DEGR- 
! EES  CELCIUS . 

0.  10086091 .25727.9. -767345.9. 78025596. 
-924  7  486589 . 6 . 98E 1 1 . -  2 . 66E 1 3 . 3 . 94E 1 4 


DATA 
DATA 
; 

READ  D< ♦ > 

! RESISTANCE  OF  PR£C 
Rd*2 . 000 
(PRINTER  IS  70 1 
SEE2 
BEE? 

INPUT 
PRINT 
BEE2 
BE20 

ASSIGNEE,  le_T0  Old* 
£NT£6  "F i  i® Itnt ( * i 


[SION  RESISTOR  :  3 /. 


“ENTER  nahc 
USING  * 2 0 >: .  ■ 


OF  -  ILE ‘ .Oldf  i  lei 
“THE  Raw  ENc  DATA 


.  :eS 


ARE  FROM  r; 


,  1  OP* 


iidf . 


FOP  i-o  TO  48 
$um«0 

FOR  J*0  TO  7 
Sijn*Sum»D(Jl*Em»(  I  i  ‘J 
NEXT  J 
TO'*  3um 
NEXT  I 

I 

J* ' 

Vo  1 f  *Em* ( 48  i 
FOR  I *5C  TO  58 

Powi  v  )  *Em»  f  I  i  *  l  Vo  i  t  -Emt  ‘III  /Ro 
Dorty  l  j  i  ■  .  03*Pow i  J) 

1  BECAUSE  OF  THE  UNCERTAINTY  IN  RP 
J  *  j  *  1 
NEXT  I 

i 

Acen*'.  .92E-4 
A ie*  *  1 . 44E -4 
Ar ig* 1 . 44E-4 
Atoo*4 . 8E-5 
Ado t* 4. BE -5 
Ato'-S. 76E-4 

T  aval  1  )  *  (  T  (  0  >  «Acen»  T(  1  >*Atop*T(2>*Arig+T<3)*Ale»»T<4)  *rtbot  i  /  A  tot 
Taval 2 ) • ( T ( 5  > *ttcen-T  <6>»HtoB*T(7i»Arig*T(8>*Alet*T(9>  »rtt>o  t ) /A  tor 
Taval3i*<T(l0i »Acen*T ( 11 >  «A  t od ♦ H 1 2  •  *Ar i a*T i 1 3 • «Ale»  »T i 14' »Abo  t ) / A  to t 
Tavg<  4 ) • i T  M  5  <  *Hcen*  T ( 1 6  i  *Atoo*  T i 1 7 ) »Ar ia*T  M  3  *  *A le* ♦ T !1 9 • *HOot i / A  to t 
Tavai  5  i  *  <  T(20i»Acen*T(2t  )-»rttoo,*Tt22**Mritp+Tt23*'*Al^**T(2FJ>  »Ado  t )  /  A  tot 
Tavai 6 ) * i T ( 25 ) *Hcen» T I  26 ) »Atoo*  T i 27 ) *Ar iq* T ( 28 • *Ale*  *T ( 29 ’ »ADot ) /Atot 


s 


rf.  f«  fu  #.*  i 


soo 
600 
6 1 0 
620 
630 
640 
650 
660 
670 
630 
630 
700 
7 1 0 
720 
730 
740 
?S0 
760 
770 
730 
730 
300 
810 
520 
330 
6  4  0 
850 

660 
'  ’A 
W  f  vJ 

330 

890 

9C0 

310 

32C 

330 

940 

350 

•J60 

370 

930 

30Q 

i  000 
'010 
1020 
'030 
:  040 
1050 
'  060 
•  0  ■’  0 
no  .  In 


ll>»  »  •«»  | 


Tavot  7  >  * i T 1 20 ) *Acen»  T ( 3 1 >  »M top*T ( 32  >  *Mr i g*  r  <  33 ' *A le* ♦ T ( 34 i «Abo  1 1 / A  t o t 
7 a *g  <  6  )  *  <  1 1  2 5  i  *Mcen  T  i  JG  >  *8 too  *  7 !  37  >  »rt  r  i  g*  7 1 33  1  *8  i  a 6  *  T  i  39  '  »Moo  t '  /  A  to  t 

T aval  9 ) • ( T ( 40 i »Hcen*T( 4 1 ) «A top* T ( 42  > *Ar i a*  7 ( 43 ) *8 1  el  *  7 t 4<i  i «Mbot ) /A tot 
i 

! UNCERTAINTY  of  tc  readings 
Dr*  .2 

Dtavq*( ( S*D t ' 2 ) ".Si/5 
7  a  *  2  0  . 

'UNCERTAINTY  OF  AMBIENT  TEMPERATURE 
Dta*  1  . 

Totr-37 .686 
Dtotr-3.8324 
Rc-54.979 
Drc*23.225 
fiea*6 .  I  09 
Dreg- . 8620 

i 

IchIP  BACK  SURFACE  TEMPERATURES 
FOR  TC  9 

T«ij)-<Pofc'ji/.083)*T((j-'  )  *5  l 

Tssum-'ssun  ji 

D  ts  .  J )  * '  1 45 . '  59-*Ddow  '  u  1  "2*Dt  '  2  !  ’  .5 

D*3avasqr*0tsava3bf»Dt3.._  i  '2 
NEXT  j 

i 

Tsava* ’ss'-in/S 

Dtsava* ( Dtsavgsqr  .5)/9 
Ota t loss* • Tsavq-7aJ/7otr 

Dato  t  *'3  tot  loss*'  (Dtsavg^2*Dt3'2;  /  t  Tsa^o- 7a  >  2*iDtorr/Totr  ■■  '2‘  '  .5 

T1*urot;o33*i 7o?f-5*oi*7a 

Dr3i**'D':otr"2*-Dr‘;'a'2’'.5 

Dmui r * t (Da tot*1 To t  r -Rea ) i '2*<GtotiOS3*0ru:*'  : "2'  ' .5 
D 1 1  ■  i Dmu  1  r " 2*0 ta ’ 2 ' * .5 

i 

! CONDUCTION  LOSS  CALCULATION 
FOR  j-i  TC  3 
Q ioss ' J ) * ■ TS'  Ji -T1 ' /Rc 

Do  . oss  i  j  •  *  3 1 oss '  J  >  » i  i  D t. s  ;  J  •  " 2  *D t  ’  ‘  2  >  t  ■  T;i  j  i  -  "  '  '■  ’  2  *  '  D r c /  Rc  :  2  1  ’  .5 
N£„"  . 

'AVERAGE  SINK  TEMPERATURE  CALCULATION 
Mud  *  4 . 32E  *  3 
At)o<un*4 . 32E-3 
A  i  n*2 . 448E "2 

I 

Tsin**<Tl4S>*T(46)*T(47l*T(48>'/4 
Dts  i  nr  -  ( ( 4»D t " 2 ' " ,5'/4 
L-2.4E-2 

PRINT  USING  "2X  .  "  "ONET  <  M  )  TAVC--TS  NU  7.IJNC.IN  HU  RA*1.E-7 

RA""  ,10a,.'" 


iDtotr/’otri "2! 


i-7(  ■■  "2* '  Drc/  Rc  : 


‘/.tJNC.  IN  HU  RA*  1  .E- 


080 

1 

090 

* 

100 

•CALCULATION  OF  ME’  °0UER. 

NU .  RH 

"10 

FOR  j-1 

TO  9 

i:c 

One  t  (  J  ' 

•Pout  J ' -Q loss ( J) 

no 

Dane  r  * ( Doou iJ>"2*Daioss(J> 

"2'-  .5 

1  4t) 

T  •  i  i  n  (  J 

t  •  >  7  avg  i  ji’Tsinlr  1/2 

15c 

Dr • iim- 

i  i  D  t  3  i  nt  "  2  *D  t  avg  "  2  i 

" .S*/2 

160 

D» 1 r • j  l 

•7 ivaijl-TsinK 

1  70 

Doe 1 1  *2 

•  Dt *  i  i« 

•* « -+  f 


■ .  ,  i  '.'i 


« 


w 

kv 


I 


I 


i 

I 


1180 
H90 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1679E- 
1200 
'  H«l, 

1310 
1320 
1  220 
1 

1  350 
1380 
•370 
1330 
’390 
1*00 
i  a  i  0 

1420 
14  30 
1  440 

1450 

1  480 
14  70 

J*r ray 

1480 
I  490 
1500 


H i J ) *One  1 1  J  i  /  t  A  to  t  »De  1 1  >  J 1 1 

Dn*H i  J  )  *  (  (  Dorn?  t  /Cine  1 1  J  >  :  '  2*  (  Doe  i  t  /De  1 1 1  J  )  )  '  22  1  '  .  5 
K ( J> * .65- 7. S94  74E -4»  f  t i Imt  J> 

Dk  *7 . 8947E -5»D t  * i lm 
K< J)*<( JJ/10. 

8no  I  J  >  *  1  .  825  -  .  00248*7 1"  i  1  m  (  J ) 

D r no* 1000 .». 0024b *[)tf  i  im 
Rho(J'*8hoiJi»'0U0. 

Cot  J)  *  .2411 !  1  *»  3. 7037E-4  )*Tf  i  ii»(  J) 

OcD*4l80.*3.7037£-4*0tt  ilm 
Cp(JI’Cpi J)*4I80. 

NIJ1-I  .  4074-2. 984E-2*rf  i  lm(J)  *3.801 8E-4*Tfi  UU)  -2-2. 7308E-6*Tfi 
-9*Tf i ln( J) *4 

Dn*D  1 1  i  i  b>*  >  8 . 735E  -  4  *  i  7 . 604E  -  4*  T(  lirBiJ)  1  *  <  2 . 4S8E -6*  3  f  :lmi  J)  ’2)  "2* 
J )  ‘  3  <  ‘  2  ’ ' .  5 
Dn*0n* 1  .  E-6 
N< J)*Nl Jl«< 1 ,£-6> 

Be  ra i J 1  * . 00248 ' 1 1 .825  - . 0024b* Tf  t lm  j  >  i 
A  1 »  a  i  J  1  *K  !  J  i  /  '  Sho  ■■  ■:  i  *03  •  j  1  > 

Da i t  a *4 ! *  a i J i ♦ < ■ Dk /K i j i ) ‘ 2* 1 Drno/Bno i ji i ' 2*  < Dcd/Cd i J  > i ' 2  > ' . 5 

^ri  ji**n  jp/fiii  a  i  j  i 


lml J >  *3*9 

1  2 . 45E -8* 


Dpr*snji»'  (  On  /  Ni  j  >  •  ‘  2*  <  0a  1 1  a/A  1  f  a  t  J )  >  '  2  1  '.5 
Nu  i  J  >  *H  (  J  I  /K  j  ' 

Dnu*Nu- j  •  » ' i Qm / A : J i i "2  *  Ok /K ( J I ) '2' ' .5 
Pernu*  '  Onu/hij .  -  •  ’  *  *  0 C  . 

C-riJl*9.3i*Beta'J>*!L'31*DeitiJ'/N(^)'2 

Dor  *0r  i  J  )  *  (  (  Dae  >  t/De  1 1 !  J  •  !  '2*<2*0r,/N(j)  i  ‘  2 )  '  .5 

Pay  i  J  i  *Gr  (  J  >  «Rf  J  > 

Dray-Ray'J)«i  t  Dgr/Gr  i  J )  i  ’2*  t  Oor/Pr  (  J )  i  '2 )  '.5 
Perray*<Dray/Rav  3)  ••100. 

Pay.2<*Say. J  /'  .5*7 

PR  INT  USING  '‘2X.2D.DDD.5'bX.2D.DD.  ;  .  /  “  I'Jne  1 1  J  i  .  De  i  1 1 ,,  i  .NijiJ)  ,Pe 
NEXT  J 

'ASSIGN  4*  ile  TO  * 

END 


iTa  T.  r.  r,- 


i  «  «»'«  t  «  W.  •  -  •  .  n  I 


LIST  OF  REFERENCES 


Chu,  R.C.,  "Heat  Transfer  in  Electronic  Systems ",  Hemisphere  Publisher  Inc., 
Heat  Transfer  1986,  v.6,  pp.  293-306. 

Knock,  R.H.,  "Flow  Visualiation  Study  of  Natural  Convection  from  A  Heated 
Protrusion  in  a  Liquid  Filled  Enclosure",  Master's  Thesis,  Naval  Postgraduate 
School,  Dec.  1983. 

Kohara,  K.  Nakao,  S.  Tsutsumi,  K.  Shibata,  H.  and  Nakata,  H.,  "High  Thermal 
Conduction  Package  Technology  for  Flip  Chip  Devices",  IEEE  Transactions  on 
Components,  Hybrids  and  Manufacturing  Technology,  v.6,  no. 3,  Sep.  1983. 

Kishimoto,  T.  Sasaki,  E.  and  Moriya,  K.,  "Gas  Cooling  Enhancement  Technology 
for  Integrated  Circuit  Chips",  IEEE  Transactions  on  Components,  Hybrids  and 
Manufacturing  Technology,  v.7,  no.3,  Sep.  1984. 

Kraus,  A.  D.,  "Thermal  Analysis  and  Control  of  Electronic  Equipment" , 
Hemisphere  Publisher  Inc.,  pp.377-407,  1983. 

Bar-Cohen,  A.  "Thermal  Design  of  Immersion  Cooling  Modules  for  Electronic 
Components",  Heat  Transfer  Engineering,  v.4,  nos.  3-4,  Jul.-Dec.  1983- 

Liu,  K.V.,  Yang.K.T.,  and  Kelleher.M.D.,  "Three  Dimensional  Natural  Convection 
Cooling  of  an  Array  of  Heated  Protrusions  in  an  Enclosure  Filled  with  a  Dielectric 
Fluid"  Proceedings  of  International  Symposium  on  Cooling  Tecnhology  of 
Electronic  Equipment,  March  17-21  1987,  pp.  486-497,  1987, 

Acharya,  S.  and  Goldstein,  R.J.,  "Natural  Convection  in  an  Externally  Heated 
Vertical  or  Inclined  Square  Box  Containing  Internal  Energy  Sources",  Journal  of 
Heat  Transfer,  November  16.  1983. 

Torok,  D.,  “Augmenting  Experimental  Methods  for  Flow  Visualisation  and  Thermal 
Perfomance  Prediction  in  Electronic  Packaging  Using  Finite  Elements",  22nd 
national  heat  transfer  conference  and  exhibition,  August  5-8  1984,  ASME  HTD 
v.32,  pp.  49-57,  1984, 

Kuhn,  D.  and  Oosthuizen,  P.H.,  "Three  Dimensional  Transient  Natural  Convective 
Flow  in  a  Rectangular  Enclosure  With  Localized  Heating" ,  Winter  Annual  Meeting 
of  ASME,  12-17  December  1986,  ASME  HTD  v.63,  pp.  55-62.  1986, 


11.  Shakerin,  S.,  Bohn,  M.  and  Loehrke,  R.I.,  "Natural  Convection  in  an  Enclosure 
with  Discrete  Roughness  Elements  on  a  Vertical  Heated  Wall" ,  Proceedings  of 
International  Heat  Transfer  Conference,  17-22  August  19S6,  v.4,  pp.  1519-1525, 
1986, 

12.  Liu,  K.V.,  Yang,  K.T.,  Wu,  Y.W.  and  Kelleher,  M.D.,  "Local  Oscillatory  Surface 
Temperature  Responces  in  Immersion  Cooling  of  a  Chip  Array  by  Natural 
Convection  in  an  Enclosure"  Proceedings  of  the  Symposium  on  Heat  and  Mass 
Transfer,  October  1-2  1987,  pp.  309-330,  1987 

13.  Park,  K.-A.  and  Bergles,  A.E.,  "Natural  Convection  Heat  Transfer  Characteristics 
of  Simulated  Microelectronic  Chips',  Journal  of  Heat  Transfer,  v.  109,  pp.  90-96, 
Feb.  1987.. 

14.  Filis,  P.  and  Poulikakos,  D.,  "An  Experimental  Study  of  the  Effect  of  Wall 
Nonuniformity  on  Natural  Convection  in  an  Enclosure  Heated  from  the.  Side",  Int. 
J.  Heat  and  Fluid  Flow,  v.7,  no.  4,  pp.  258-265,  Dec.  1986, 

15.  Anderson,  R.  and  Bohn,  M.,  "Heat  Transfer  Enhancement  in  Natural  Convection 
Enclosure  Flow",  Winter  Annual  Meeting  of  ASME,  9-14  December  1984  Journal 
of  Heat  Transfer,  v.lOS,  no.  2,  pp. 330-336,  May  1986. 

16.  Yang,  R.  and  Yao.L.S.,  " Natural  Convection  Along  a  Finite  Vertical  Plate", 
Journal  of  Heat  Transfer,  v.109,  pp.  363-413,  May  1987, 

17.  Beckermann,  C,  Ramadhyani,  S.  and  Viskanta,  R,  “Natural  Convection  Flow  and 
Heat  Transfer  Between  a  Fluid  Layer  and  a  Porous  Layer  Inside  a  Rectangular 
Enclosure",  Journal  of  Heat  Transfer,  v.I09,  pp.  363-413,  May  1987, 

18.  Kulkarni,  A.K.,  Jacobs,  H.R.  and  Hwang,  J.J.,  " Similarity  Solution  for  Natural 
Convection  Flow  over  an  Isothermal  Vertical  W all  Immersed  in  Thermally  Stratified 
medium",  Int.  J.  Heat  Mass  Transfer,  v.30,  no.  4,  pp.  691-698,  1987, 

19.  Lankford,  K.E.  and  Bejan,  A.,  "Natural  Convection  in  a  Vertical  Enclosure  Filled 
With  Water  near  4  °C\  Journal  of  Heat  Transfer,  v.108,  pp.  755-763,  Nov.  1986 

20.  Bohn,  M.S.  and  Kirkpatrick,  A.T.,  "High  Rayleigh  Number  Natural  Convection  in 
an  Enclosure  Heated  from  Below  and  From  the  Sides",  Natural  Heat  Transfer 
Conference,  24-28  July  1983  ASME  HTD,  v.26,  pp.  27-33,  1983. 

21.  Kline,  S.J. and  McClintock,  F.A.,  "Describing  Uncertainties  in  Single-sample 
Experiments" ,  Mechanical  Engineering  ,  p.3,  Jan.  1953. 

22.  Holman,  J.P.,  "Heat  Transfer",  Mc-Graw  Hill,  pp.  265-305,  1981. 


>VI  *  * 


r 


.*§  »>*' i 


23.  Kays..  W.M.  and  Crawford.,  M.E.,  Convective  Heat  and  Mass  Transfer"  .  Mc- 
Graw  Hill,  pp.  313-331,  1980- 


i 

f 

i 


w 


m 


.‘i*. 


*  u* •  i,#  Mji’ji*  y.y 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  VA  22304-6145 

Library.  Code  0142 
Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

Professor  \1.  D.  Kelleher,  Code  69KK 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Y.  Joshi,  Code  69JI 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Department  Chairman  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Deniz  Kuvvetleri  K.ligi 
Personel  ve  Egitim  Daire  Baskanligi 
Bakanliklar;  AN  K  A  RA-T  urkey 

Istanbul  Technical  University 
Makina  Muhendisligi  Fakultesi 
Isi  Transferi  Kursusu 
.Ylacka/ISTANBUL-Turkey 

Bosphorus  University 
Makina  Muhendisligi  Fakultesi 
Isi  Transferi  Kursusu 
I  ST  ANBU  L-T urkey 

Middle  East  Technical  University 
Makina  Muhendisligi  Fakultesi 
(si  Transferi  Kursusu 
ANKARA-Turkey 

Turgay  Pamuk 
Ogretmenler  Sitesi 
Odul  Sok.,  A  Blok,  No  -  II 
Ambarli/ISTANBUL-Turkey 


| 

I 

s 


|| 

& 


